TAILS OF THE ENDPOINT DISTRIBUTION OF DIRECTED 

POLYMERS 

JEREMY QUASTEL AND DANIEL REMENIK 



Abstract. We prove that the random variable T = argmax teH {.42(£) — t 2 }, where Ai 
is the Airy2 process, has tails which decay like e~ ct . The distribution of T is a universal 
distribution which governs the rescaled endpoint of directed polymers in 1 + 1 dimensions 
CN| for large time or temperature. 

o 

_ 1. Introduction 

Qh 

Consider the following model of a directed polymer in a random environment. A polymer 

path is a nearest-neighbor random walk path it = (ttq,tti, . . . ) in Z started at the origin, 

that is, 7To = and iTk — ftk-i = ±1- On Z + x Z we place a collection of independent 

random weights {<^i,j} >(] - Gj - The weight of a polymer path segment it of length N is 

defined as 
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for some fixed /3 > which is known as the inverse temperature. If we restrict our attention 
to paths of length TV which go from the origin to some given x G Z then we talk about a 
l — ' point-to-point polymer, defined through the path measure 

Q^ for 7r of length AT going from the origin to x and Q^ 1 ^ (it) = otherwise. The normalizing 

constant Z pomt (N,x) = X^7r-7r(o)=o n(N)=x ^Jv(tt) is known as the point-to-point partition 
function. Similarly, if we consider all possible paths of length iV then we talk about a 

^j point-to-line polymer, defined through the path measure 



QN e W = ^V)^(TT) 



for 7r of length N and Qjf c (n) = otherwise, with the point-to-line partition function 
Z linc (N) = J2k=-N Z point (N, k). 

Our main interest will be the point-to-line case. The main quantities of interest in this 
case are the partition function and the position of the endpoint of the randomly chosen 
path, which we will denote by T/v- It is widely believed that these quantities should satisfy 
the scalings 

log(Z linc (A0) ~ aN + WV 1/3 x, and T N ~ N 2 ' ? 'T, 

where the constants a and b may depend on the distribution of the tOij and /3, but x an d 
T should be universal (up to some moment assumptions on the cjjj's). 

While there are few results available in the general case described above, the zero- 
temperature limit (3 — > oo, known as last passage percolation, is very well understood, 
at least for some specific choices of the environment variables ojij. We will restrict the 
discussion to geometric last passage percolation, where one considers a family {coij}i^z+,jez 

l 
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of independent geometric random variables with parameter q (i.e. P(a^ ,• = k) = q(l — q) k 
for k > 0) and defines the point-to-point last passage time by 

N 

L(N, y) = max > w» «.(,-\ 

and the point-to-line last passage time by 

L(N)= max L(iV,y). 

y=—N,...,N 

We remark that this model is usually defined on (Z + ) 2 , which corresponds to rotating our 
picture by 45 degrees and working on the dual lattice. Although the exact results we will 
describe next have been proved for that case, the picture in our situation is morally the 
same, and hence for simplicity we present the results for last passage percolation on Z + x Z. 

We define the rescaled process t h-> H^{t) by linearly interpolating the values given by 
scaling L(N, y) through the relation 

L(N, y) = Cl N + c 2 N 1 / 3 H N (c 3 N- 2 / 3 y), 

where the constants q have explicit expressions which depend only on q and can be found 
in [Joh03]. The point-to- line rescaled process is then given by 

G(N) = sup H N (t), 

tel-ca-ATVVsiVVS] 
and it is known in this case [BR01] that 

(1.1) G(N) ~ aN + bN^ 3 x 

with x having the Tracy- Widom largest eigenvalue distribution for the Gaussian Orthogo- 
nal Ensemble (GOE) from random matrix theory [TW96] (the analogous result holds in the 
point-to-point case with x now having the Tracy- Widom largest eigenvalue distribution for 
the Gaussian Unitary Ensemble (GUE) [TW94]). On the other hand, Johansson [Joh03] 
showed that 

H N (t) ->A 2 (t) -t 2 

in distribution as N — > oo, in the topology of uniform convergence on compact sets. Here 
A 2 is the Airy2 process, which we describe below, and which is a universal limiting spa- 
tial fluctuation process in such models. As a consequence of Johansson's result (see also 
[CQR12]), (1.1) translates into 

(1.2) Pfsup (A 2 (t) - t 2 ) <m)= F GOE (4 1/3 m) 

(the 4 1 ' 3 arises from scaling considerations, or alternatively from the direct proof given in 
[CQR12]). 

The Airy2 process was introduced by Prahofer and Spohn [PS02], and is defined through 
its finite-dimensional distributions, which are given by a Fredholm determinant formula: 
given xo, . . . , x n € M. and to < ' ' ' < t n in R, 

(1.3) F{A 2 (t ) < x , . . .,A 2 {t n ) < x n ) = det(I - f 1/2 Kcxtf 1/2 )L2 ({io ,..., in }xR), 

where we have counting measure on {to, ■ • • ,t n } and Lebesgue measure on M, f is defined 
on {to, . . . ,t n } x M by f(tj,x) = l X £(x -,oo), an d the extended Airy kernel [PS02; FNH99; 
Mac94] is defined by 

K Hem J/o^rfAe-^-^Ai^ + ^Ai^ + A), if t > t> 

extl '^' ' U "\-/_° oo( iAe- A ('-'')Ai(e + A)Ai(e' + A), if t < f, 

where Ai(-) is the Airy function. In particular, the one point distribution of A 2 is given by 
the Tracy- Widom GUE distribution. An alternative formula for A 2 due to [PS02], which 
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is the starting of the proofs given in [CQR12] of (1.2) and (1.8) below, and also of the main 
result of this paper, is given by 

(1.4) F(A 2 (to)<xo,...,A 2 (tn)<x n ) 

= det(l - K Ai + P Xo e^-^ H P Xl e^-^ H ■ ■ ■ P Xn e^-^ H K M 

where K^ is the Airy kernel 

,o 

K Ai (x, y) = dX Ai(x - A) Ai(y - A), 

J — oo 

H is the Airy Hamiltonian H = —d 2 + x and P a denotes the projection onto the inter- 
val (— oo, a]. Here, and in everything that follows, the determinant means the Fredholm 
determinant on the Hilbert space L 2 (M), unless a different Hilbert space is indicated in 
the subscript (the last formula (1.4) should be compared with (1.3), where the Fredholm 
determinant is computed in an extended space). The equivalence of (1.3) and (1.4) was 
derived in [PS02] and [PS11]. We refer the reader to [CQR12; QR12a] for more details. 

Coming back to geometric last passage percolation, we turn to the random variables 

T/v = inf \t: supH^is) = sup-ffjv(s)}, 

8<t sSK 

which correspond to the location of the endpoint of the maximizing path with unconstrained 
endpoint (that is, the zero-temperature point-to-line polymer). From the above discussion 
one expects the following: 

Theorem 1. Let T = argmax ieK {„42(£) — t 2 }. Then, as N — )■ oo, T/v -^ T in distribution. 

This result was proved by Johansson [Joh03] under the additional hypothesis that the 
supremum of A%(t) — t 2 is attained at a unique point. The uniqueness was proved, using 
two different methods, by Corwin and Hammond [CH11] and by Moreno Flores and us 
[MQR12]. 

Although the result of Theorem 1 has only been proved in the case of geometric (or ex- 
ponential) last passage percolation, the key point is that the polymer endpoint distribution 
is expected to be universal for directed polymers in random environments in 1 + 1 dimen- 
sions, and even more broadly in the KPZ universality class, for example in particle models 
such as asymmetric attractive interacting particle systems (e.g. the asymmetric exclusion 
process), where second class particles play the role of polymer paths. This problem has 
received quite a bit of recent interest in the physics literature, see [MQR12] and references 
therein for more details. 

In [MQR12] we obtained an explicit expression for the distribution of T ■ More precisely, 
we obtained an explicit expression for the joint density of 

T = argmax{„4 2 (£) - t 2 } and M = max{„4 2 (i) - t 2 }, 

which we will denote as f(t,m). To state the formula we need some definitions. Let B m 
be the integral operator with kernel 

B m (x,y) = Ai(x + y + m). 

Recall that Ferrari and Spohn [FS05] showed that Fgoe can be expressed as the determi- 
nant 

(1.5) F GOE (m) = det(I - P B m P ), 

where P a denotes the projection onto the interval [a, oo) (the formula essentially goes back 
to [Sas05]). In particular, note that since i 7 GOE( ?Tl ) > for all m € R, (1.5) implies that 
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I — PoB m Po is invertible. For t,m£R define the function 

^ t , m (x) = 2e xt [t Ai(x + m + t 2 ) + Ai'(x + m + t 2 )] 

and the kernel 

*t, m (x,y) = 2 1 / 3 V i , m (2 1 / 3 x)V- i , m (2 1 / 3 y). 
Then the joint density of T and M is given by 

/(*, m) = det(J - P B 4l/3m P + P^ t ,mPo) ~ F GOE (4 1/3 m) 

= tT[(I-P B±y 3m P )- 1 P y t>m P ]F GO v(4 1/3 m). 

Integrating over m one obtains a formula for the probability density f e nd(t) of T, although 
it does not appear that the resulting integral can be computed explicitly. One can readily 
check nevertheless that f e nd(t) is symmetric in t. Figure 1, taken from [MQR12], shows a 
plot of the marginal T density. 

The goal of this paper is to study the decay of the tails of T. We will prove: 
Theorem 2. There is a c > such that for every k > ^- and large enough t, 

e -^ < P (|r| > t ) < ce -§t3 +24 2 +0(t 3/ 2) ^ 

We believe that the correct exponent is the — , obtained in the upper bound (we remark 
that we have not attempted to get better estimates on the lower order terms in the upper 
bound). The tail decay of order e~ ct confirms a prediction made in the physics literature 
in Halpin-Healy and Zhang [HHZ95], see also Mezard and Parisi [MP92]. Their idea is to 
argue by analogy with the argmax of Brownian motion minus a parabola. In that case one 
has a complete analytical solution [Gro89]. 

We will give two proofs of the upper bound, both in Section 2. The first one is based on a 
direct application of the formula (1.6) for the joint density of T and Ai. The second proof 
will start from a probabilistic argument and then use the continuum statistics formula for 
the Airy2 process obtained in [CQR12] to estimate the probability that the maximum is 
attained very far from the origin. This formula corresponds to the continuum limit of (1.4) 
and is given as follows (see [CQR12] for more details). Fix £ < r and g £ H l ([l, r\) and 
define an operator @? ir i acting on L 2 (M) by ©fg r i/(') = u(r, •), where u(r, •) is the solution 
at time r of the boundary value problem 

d t u + Hu = for x <g{t), t £ (£,r) 

(1.7) u(£,x) = f(x)l x<g(e) 

u(t, x) = for x > g{t). 

Then 

(1.8) ¥(A 2 (t) < g(t) for t G [£,r]) = det(/ - K Ai + ef^e^^Ai) • 

We remark that in the second proof actually get an upper bound with a larger 0(t 2 ) 
correction in the exponent. 

Not surprisingly, the lower bound turns out to be more difficult (in fact, for the upper 
bound we can basically use the estimate | det(/ + ^4) — det(/ + i?)| < \\A — B\\ie 1+ " A " 1+ " B " 1 
for trace class operators A and B directly to estimate the decay by computing the trace 
norm of two operators; no such estimate is available for the lower bound). In this case we 
will have to use a probabilistic argument to extract the lower bound from the well-known 
exact asymptotics for the tails of the GUE distribution, and then show that the remaining 
terms are of lower order. For this last task we will use again (1.8), but the argument is 
much more complicated than for the upper bound. Interestingly, it will involve turning 
an instance of (1.8) which mixes continuum and discrete statistics for Ai back into an 
extended kernel formula. 
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Remark 1.1. 

1. A few days before submitting this article, we became aware of the very recent work 
of Schehr [Schl2], where he obtains, using non-rigorous arguments, an alternative 
formula for the joint distribution function of A4 and T ■ His formula is obtained by 
taking the limit in N of a known formula for the joint distribution of the maximum 
and location of the maximum for the top line of N non-intersecting Brownian ex- 
cursions, which is expected to converge to the Airy2 process. The resulting formula 
is expressed in terms of quantities associated to the Hastings-McLeod solution of 



the Painleve II equation, and has tails decaying like e 



2. During the refereeing process, Baik, Liechty, and Schehr [BLS12] proved the equiva- 
lence of the formula of [Schl2] and (1.6). Hence the rigorous validity of the formula 
of [Schl2] is established based on [MQR12], as well as the tail decay. 

Acknowledgements. The authors would like to thank the referee for a careful reading of 
the article. Both authors were supported by the Natural Science and Engineering Research 
Council of Canada, and DR was supported by a Fields-Ontario Postdoctoral Fellowship. 

2. Upper bound 

Throughout the paper c and C will denote positive constants whose values may change 
from line to line. We will denote by || • || op , || • ||i and || • ||2 respectively the operator, trace 
class and Hilbert-Schmidt norms of operators on L 2 (M) (see Section 3 of [CQR12] for the 
definitions or [Sim05] for a complete treatment). We will use the following facts repeatedly 
(they can all be found in [Sim05]): if A and B are bounded linear operators on L 2 (M), then 



(2.1) 



\AB\\i < \\A\\ 2 \\B\\2, \\AB\\ 2 < ||A|| p||5|| 2 , ||AB|| 2 < P|| 2 ||.B|| p, 

Pllop < ||A|| 2 < P||i, 



\A\\l = / dxdyA(x,yf 

JR2 



where in the last one we are assuming that A has integral kernel A(x,y). We will also use 
the bound 

(2.2) |det(/ + A) - det(I + B)\<\\A- B\\ 1 e^ 1+ ^ B ^ +1 < \\A - B|| ie P-B||i+2||.B||i+i 

for any two trace class operators A and B. 

We recall that the shifted Airy functions <j)\{x) = Ai(x — A) are the generalized eigenfunc- 
tions of the Airy Hamiltonian, as H<p\ = X(p\, and the Airy kernel i^Ai is the projection 
of H onto its negative generalized eigenspace (see Remark 1.1 of [CQR12]). This implies 
that e s Kai has integral kernel 

poo 

(2.3) e sH K Ai (x,y) = d\e~ sX Ai(x + A) Ai(y + A). 

Jo 

It also implies that e aH K^ie bH Ka\ = e( a+b > H Kx\- We will use this fact several times in 
this and the next section. 

2.1. First proof. We start by writing 

P(T >t)< P(T > t, M > -2t) + F(M < -2t). 

By (1.2) the second probability on the right side equals -Fgoe( — 2 5 > 3 t) < ce~3* , where the 
tail bound can be found in [BBD08]. Thus it will be enough to prove that 



(2.4) P(T > t, M > -2t) < ce 



.4/3 
3 



t3 +2t 2 +C , (i 3/ 2) 
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Let s > t. We will assume for the rest of the proof that m > —2b. Using (1.5) and the 
first formula in (1.6) we get from (2.2) that 
(2.5) f(s,m) < ||P ^ irnJ Po||ie 1+2||PoB 4V3 m fl»||i+||Po^, m Po||i > 

Using the identity 

/oo 
du Ai(a + u)Ai(6 - u) = 2~ 1/3 Ai(2~ 1/3 (a + &)) 
-oo 

and letting e = i _1 we may write 

P B 4l/3m P = 2 1 / 3 Q 1 Q 2 with Qx(x, A) = l x > Ai(2 1 / 3 x + m + A)e5 e \ 



(2.6) 



-he\ 



Q 2 (A ; y) = e-3 £A Ai(2 i / d y + m - A)l v > . 



V3, 



Lemma 4.2 now gives 

(2.7) ||P S 4 i/3 m Po||i < ci 3 / 2 . 

On the other hand, recall that the trace norm of an operator ^ acting on L 2 (M) is defined 



as 



W\ 



X)( e "'i*i 



n=l 



where {e n } n >i is any orthonormal basis of L 2 (M) and |vl/| = vW is the unique positive 
square root of the operator *&**$/. For the case ^f = Pq^s^Pq, since ^ is a rank one 
operator it is easy to check that \I/*\I/ has only one eigenvector, and in fact it is given by 
lx>o^- s ,m(2 1/3 x) with associated eigenvalue A s , m = 2V 3 ||P ^,m(2 1/3 OII|||^oV'- s ,m(2 1/3 -)||i. 
We deduce that H-Po^s.m-Fblli = v/A s>m , and then by (2.5) and (2.7) we get 



(2. 



-it 



dmf(s,m)< dm^X^e 1+ct3/2+ ^ 

J-2t 



Now Lemma 4.3 gives 



dm y%^ = 2 1 / 6 / dm ||P V's, m (2 1/3 -)||2||PoV'-s,m(2 1/3 -)||2 

2t J-2t 

1/2 



<2 l/6 



dm||P V's,m(2 1/3 -) 



2t 

< ce-l s3+2st > 



2/ 



dm||P V'- s , m (2 1/3 -) 



nl/2 



and it is not hard to see from the proof of Lemma 4.3 that A Sjm is bounded uniformly for 
m > —2t, s > t and large enough t. We deduce then from (2.8) that J_ 2t dm f(s,m) < 



ce 3 s + 2si +C(t ^ anc j h ence 

/'OO /'OO /'OO 

F(T > t, M > -2t) = ds I dmf(s,m)<c dse'z 

Jt J-2t Jt 



s 3 +2st+<D(t 2 ) 



< ce s 



4 t 3 +2i 2 +C , (t 3/ 2) 



where the last estimate can be easily obtained from an application of Laplace's method, 
see the proof of Lemma 4.1 for a similar estimate. This gives (2.4) and the upper bound 
of Theorem 2. 

2.2. Second proof. Since we already have a full proof of the upper bound, we will skip 
some details. The key result for this proof is the following 

Proposition 2.1. Fix L > 1. Then there is a c > such that for every m > 



sup A 2 {x) > m + 1 < ce ' ; 

,xe[-L,L] J 



i m .3/2 
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Figure 1. Plot of the density of T compared with a Gaussian density with 
the same variance 0.2409 (dashed line). The excess kurtosis E(7~ 4 )/IE(T~ 2 ) 2 — 
3 is -0.2374. 



Proof. By (1.8) we have, writing g m (s) = s 2 + m, 



sup A 2 (x) < m + 1 J > PI sup (A 2 (x) - x 2 ) < m) 

\xe[-L,L] J \xe[-L,L] J 

= det(/ - K Ai + e LH K Ai e^ L L f LH K M ) , 



(2.9) 



where we have used the cyclic property of the determinant together with the identity 
e 2LH Km = e LH Kmz LH Km (see the remark after (2.3)). Now recall from Theorem 1.3 of 
[CQR12] that 

F GOE (4 1/3 m) = det(J - e LH K Ai R L e LH Km) , 

where Rl is defined in (1.5) in [CQR12]. Therefore by (2.2) we deduce that 
(2.10) 
det(/ - Km + e LH K M ep LM e LH KM) > F GOE (^m) - \\A - S || ie i+ll^-B|| 1 +2||B|| lj 



where 



A = Km ~ e LH K M Q 9 r L L] e LH K M and B = e LH K A ,R L e LH K M 



and we have used the triangle inequality in the exponent. Now ||i?||i can easily be bounded 
by some constant uniformly in m > by an argument similar to the one used to obtain 
(2.7). On the other hand, using the decomposition of 8j_V ri given in (3.4) in [CQR12] we 
have 

A - B = e LH K Ai tt L e LH K AU 

where Vt L = (R L - P m+L 2R L P m+L 2) - (e~ 2LH - P m+L 2e- 2LH P m+L2 ). By Lemma 4.4 we 
get 11^4 — B\\i < ce~ vm for some rj > ,, and then using this and (2.10) in (2.9) we obtain 



sup A 2 {x) < m + 1 > F GOE (4 1/3 m) -ce""™ 3 ' 2 . 
v *e[-i,i] / 



The result follows from this and the asymptotics [BBD08] i ? GOE(4 1 ' 3 m) > 1- 



-3/2 ~^m 3 / 2 

-cm i e 3 

□ 
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Using Proposition 2.1 we can derive the upper bound. Start by observing that, for fixed 
ae(0,l), 

(2.11) 

P(Te [8,8 + 2]) <P( sup {A 2 (x)-x 2 ) >A 2 (0)\ <F\ sup A 2 (x) > A 2 {0) + s 2 ) 
\xe[s,s+ 2 ] J \xe[s,s+2] J 

<P( sup A 2 (x)> (l-a)s 2 ) +P(A 2 (0) <-as 2 ) . 

\x£[s,s+ 2 ] J 

The last probability equals Fq\je(—cts 2 ), and then the asymptotics obtained in [BBD08] 
give PL42(0) < —crs 2 ) < ce~T2°" s . For the other term on the right side of (2.11) we use 
Proposition 2.1 and the stationarity of the Airy2 process to write, for s > t, 

p( sup A 2 (x)>(l-a)s 2 ) < ce -f[d-^ 2 -i] 3/2 . 

Therefore 

P(Te [8,a + 2]) < ce -™n{^ 3 s 6 4l(i-<?)s 2 -i] 3/2 }. 

Take <r = 4 2 ' 3 s so that the minimum in the above exponent equals n[s 2 — 4 2 ' 3 s — l] 3 ' 2 = 
|s 3 + 0(s 2 ) for large enough s. We deduce that 

¥(Te[s,s + 2])<ce-l s3+ °^. 
Summing this inequality over intervals of the form [t + 2k, t + 2(k + 1)] for k > gives 

¥{T>t)<ce-l t3+ °^ 

for some c > and large enough t, and now the upper bound in Theorem 2 (with a worse 
0(t 2 ) correction) follows from the symmetry of T. 

3. Lower bound 

As we mentioned, the lower bound turns out to be more delicate, because in this case 
a simple bound like (2.2) is not available. The main idea of the proof is to compare the 
probability we are interested in with an expression involving the one-dimensional marginal 
of the Airy2 process, and then extract the lower bound from known asymptotics for the 
Tracy- Widom GUE distribution. This comparison will introduce an error term, and most 
of the work in the proof will be to show that this error term is of lower order. 

The first step in the comparison is to write, for t > 0, /3 > and s > 0, 

(3.1) P(|T| >t)> F(A 2 (x) - x 2 < f3t 2 Vx < t, A 2 {t + s)-(t + s) 2 > fit 2 ) . 

The idea is the following. If s is now taken to be reasonably large, then A 2 (t+s) and A 2 (x), 
x < t, should have decorrelated somewhat. Assuming they have completely decorrelated, 
we would have that the right side is bounded below by P(A 2 (t + s) — (t + s) 2 > fit 2 ) , which 
has the correct decay if we choose s = at. So the whole proof comes down to estimating 
the correction coming from the correlation. 

Of course, from (3.1) we have 

(3.2) P(|T| > t) = F(A 2 (x) - x 2 < fit 2 Vi < t) 

- F(A 2 (x) -x 2 <pt 2 \/x<t, A 2 (t + s)-(t + s) 2 < f3t 2 ) , 
The bound will then follow from the following 
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Lemma 3.1. Let (3 > 3. There is an ao > (which depends on (3) such that if a £ (0, ao) 
and s = at, then for large enough t we have 

F(A 2 (x) - x 2 < (3t 2 Vx < t, A 2 (t + s)-(t + s) 2 < (3t 2 ) 

< F(A 2 (x) - x 2 < f3t 2 V x < t) F(A 2 {t + s)-{t + s) 2 < (3t 2 ) 

.[l + l a2 ( i + s )-3 e -!(/3+D 3 / 2 (^) 3 ] ; 

where a 2 is defined implicitly in (3.3). 

To see how the lower bound follows from this, let f3 > 3 and choose a as in the lemma. 
Then letting s = at and using the lemma and (3.2) we get 

P(|T| > t) > F(A 2 (x) - x 2 < pt 2 Vx < t) 

■ (l - F(A 2 (t + s)-(t + s) 2 < f3t 2 ) [1 + \a 2 (t + s )-3 e -|(^+i) 3 / 2 (t+ s )^ . 

Now if we let p = F(A 2 (x) - x 2 < Vx G R) = F GO e(0) > then, since t > 0, the first 

probability on the right side above is larger than p$. On the other hand, 

(3.3) 

F(A 2 (t + s)-(t + s) 2 < pt 2 ) = F GUE ((/3 + l)t 2 + 2st + s 2 ) < F GUE ((/3 + l)(i + s) 2 ) 

<l_ a2 ( t + s) -3 e -!(/3+i) 3 / 2 ( i+ s) 3 

for some explicit constant a 2 > and large enough t, see for instance [BBD08] for the 
precise bounds on the tails of the GUE distribution. This implies that 

F(\T\ >t)> Po(^a 2 (t + s )-3 e -i(/m) 3/2 (*+*) 3 + l ffl 2 (f + s) -6 e -|(/3 + l)3/ 2 ( t+s )3^ 

>c(l + a)- 3 t- 3 e-f( 1+ ^ 3/2 ( 1+a ) 3t3 , 

and now the lower bound in Theorem 2 follows from choosing (3 > 3 and a small enough 
a > so that f (1 + /3) 3 / 2 (l + a) 3 < k. 

Our goal then is to prove Lemma 3.1. For this we need an expression for the probability 
we want to bound. The answer follows from a simple extension of the result in [QR12b], 
where we obtained an explicit expression for probabilities of the form 



sup(^,2(x) — x ) < m 
y x<t y 

and showed that they to correspond (after a suitable shift) to the one-dimensional marginals 
of the Airy2_>.i process [BFS08]. To state the extension of that formula, define for o,t£l 
the operators 

Qa,tf(x) = f(2(a + t 2 ) - x) and M 0)t /(x) = e 2t ^- a '^ f(x) 

acting on / £ L 2 (IR). We will say that an operator acting on L 2 (M) is identity plus trace 
class if it can be written in the form I + A with A a trace class operator. 

Lemma 3.2. With the above definitions, and for any a,b EM and s > 0, 

(3.4) F(A 2 (x) - x 2 < a Vx < t, A 2 (t + s) - (t + s) 2 < b) 

= det(J - K Ai + K M (I - M a , t0ayt )P a+t 2e- sH P b+{t+s) 2e sH K Ai ) . 
Moreover, the operator inside this determinant is identity plus trace class. 

Proof. For L > it is straighforward to adapt the proof given in [CQR12] of the continuum 
statistics formula (1.8) to deduce that 

F(A 2 (x) -x 2 < a VxG [~L,t], A 2 {t + s) - (t + s) 2 < b) 

= det(7 - K Ai + e { -LA^ sHp b+{t+s?z {L+t+s)H KM 
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where ©r-z,t] is defined as ©f_ L tl (see (1.7)) for g(x) = x 2 + a. Since e^+'+^i^Ai = 

e sH K Al e {L+t ^ H K Ai and #Ai = e sIi K Ai e~ sH K Al = e~ sH K Ai e sH K Ai (see the remark after 
(2.3)), we can use the cyclic property of the determinant to turn the last determinant into 

(3.5) det(l - K M + e^ L+ ^ H K Ai Q { _ LA e- sH P b+{t+s)2 e sH K A ^j . 
Rewriting the above operator as 

(3.6) /- e - sH ^Ain + ( t+ ,) 2 e sH ^Ai-e( L+< ^K Ai (e-( L+4 ) H -G [ _ L ,^e-^n + ( i+s)2 e^^Ai, 

it follows from an easy adaptation of the proof of Proposition 3.2 of [CQR12] that the 
operator is identity plus trace class. On the other hand, in the proof of Theorem 1 of 
[QR12b] it was shown that 

e (L +i )H KA . e y Ka . { j _ M a j 6 a,t)Pa+ti 

1 ' L— >oo 

in Hilbert-Schmidt norm, and a straightforward extension of the proof shows that the same 
holds if we post-multiply both sides by e~ sH P b+ ( t+s }2N, where Nf(x) = (1 + x 2 ) 1 ' 2 f(x). 
Thus, since N~ l e sH K M is Hilbert-Schmidt by (3.3) in [CQR12] deduce by (2.1) that 

e (L + t)H KA . e[ _ Lit] e- 8H P bHt+3) 2e sH K Ai -£— ► K Ai (I-M a jg a , t )P a+t2 e- sH P bHt+s)2 e sH K Ai 

in trace class norm. This together with (3.5) and (3.6) yields (3.4) and, in particular, the 
fact that the operator inside this determinant is identity plus trace class. □ 

The key to obtain Lemma 3.1 from Lemma 3.2 is to turn the last determinant into the 
determinant of a 2 x 2 matrix kernel. Observe that, since M a j and P a +t 2 commute and 
P a +t 2 Qa,t = Qa,tP a +t 2 -> the formula (3.4) can be rewritten as 

(3.7) F(A 2 (x) - x 2 < a Vx < t, A 2 (t + s) - (t + s) 2 < b) 

= det(/ - K Ai + K Ai (I - Q)e~ sH {I - P 2 )e sH K Ai ) 

where 
(3.8) 
Pi = P a+t 2, P 2 = P b+ ( t+ s) 2 , an d Q = Pi(I + M a)t Q ajt ) = P a+t 2 + M ayt Q a jP a+t 2 . 

Note that Q 2 = Q (although Q is not a projection in L 2 (R), as it is an unbounded operator). 
This formula has exactly the same structure as the formula (1.4) for the finite-dimensional 
distributions of the Airy2 process (in the case n = 2) which, we recall, is equivalent to the 
extended kernel formula (1.3). The equivalence of the two types of formulas was developed 
by Prahofer and Spohn [PS02] and Prolhac and Spohn [PS11], and later made rigorous 
and extended to the Airyi case by us in [QR12a]. The same method will work for (3.7), 
and the result is the following: 

Proposition 3.3. With the notation introduced above, 

(3.9) F(A 2 (x) - x 2 < a Vx < t, A 2 (t + s) - (t + s) 2 < b) 

Qi^AiPi Qe- sH (K M -I)P 2 
P 2 e sH K M P 1 P 2 K Ai P 2 



det [I - r 



r- 1 

' LH 



with Gf(x) = e- 2tx cp- 1 (x)f(x) and <p(x) = (1 + x 2 ) 1/2 . 



where 

" G 
G 

In particular, the operator above is identity plus trace class. 

The conjugation by V is needed to make the operator trace class. Note that there is a 
slight difference between this formula and the ones for the Airyi and Airy2 processes: the 
formula is not written in the most symmetric way, as the first column in the brackets is 
post-multiplied by P\ instead of Q. Formally there is no difference, because since Q 2 = Q 
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one can pre-multiply the whole matrix by 



Q 
Q 



and then use the cyclic property of the 



determinant to turn this into a post-factor of P\Q = Q for the first column. But this form 
of the formula will turn out to be better for obtaining the desired bounds in Lemma 3.4. 

The proof of Proposition 3.3 follows the steps of the proofs in [PS11; QR12a], but given 
the slight difference noted above, and since it is short and easy to present in the two- 
dimensional case, we include the details. 

Proof of Proposition 3. 3. We remark that all the manipulations performed on Fredholm 
determinants below rely on knowing that the operators inside each of them are identity 
plus trace class (see the proof of Theorem 1 of [QR12a] for more details), but this can be 
seen in each case from Lemma 3.2 and Lemma 3.4 and similar estimates. 

Let D denote the determinant in (3.9). The operator inside it can be factored as 
/ GQe- sH P 2 G- 1 



I 



QK M Pi ~ Qe~ sH P 2 e sti K M Pi Qe~ sH K Ai P 2 - Qe~ sH P 2 K Ai P 2 
P 2 e sH K Ai P 1 P 2 K Ai P 2 

Note that the determinant of the first matrix on the right side above is 1 , so D equals 

-sH -D_„aH ir . . e> n„-sH xr . ,r>_ n n -sH , 



detfl 



QK M Pi - Qe~ sH P 2 e sti K Ai Pi Qe' sH K Ai P 2 - Qe~ sti P 2 K M P 2 
P 2 e sii K Al P x P 2 K Ai P 2 

QK Ai -Qe- sH P 2 e sH K Ai 
P 2 e sH K Ai 



L 2 ( 



det(/-r 



Pi e~ sH P 2 











LH 



Since K A \ is a projection we may pre-multiply each entry in the second bracket by K A \ 
and then use the cyclic property of the determinant to get 



D 



detfl 
det( I 



KaiPi K M e- sH P 2 











QK Ai -Qe- sH P 2 e sH K Ai 
P 2 e sH K Ai 



l'H 



KmQKm - K Ai Qe- sH P 2 e sH K Ai + e~ sH K M P 2 e sH K M 





LH 



= det(J - K Ai QK Al + K Ai Qe~ sH P 2 e sH K Ai - e~ sH K Ai P 2 e sH K Ai ) L2(R) . 
This last determinant equals the one on the right side of (3.7), and the result follows. □ 

The 2x2 matrix kernel formula is useful because it will allow us to extract easily the 
first two factors in the bound given in Lemma 3.1. This idea was introduced by Widom 
[Wid04], where he studied the asymptotics in t of F(A 2 (0) < si, A 2 (t) < s 2 ). 

Proof of Lemma 3.1. We start with the formula in Proposition 3.3, with a = b = j3t 2 , and 
use the idea introduced in [Wid04]: factor out the two diagonal terms in the determinant 
and then estimate the remainder. More precisely, we write 

Qe~ sH (K P 



I-T 

I - I' 



QK M Pi 

P 2 e sH K Al P x 



, Ai -I)P 2 1= ( QK M Pi 

P2K Ai P 2 \ \ [ P 2 K Ai P 2 

(I-QK Ai P l )- 1 Qe- sH (K Ai -I)P 2 

(I - P^mP^P^KmPi 



and then recognize that the determinant of the first factor on the right side equals 

det(/ - GQK Xl PiG- 1 ) det(/ - GP 2 K Al P 2 G~ l ) 

= F{A 2 {x) - x 2 < Vx < t) ¥{A 2 {t + s)-{t + s) 2 < 0) . 



det J-T 
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To get the last equality, observe first that det(I — GP 2 K A iP 2 G~ 1 ) = -Fgue((£ + s) 2 ) which 
is the second factor on the right side. For the first one we note that, by the cyclic property 
of determinants and the facts that K Ai = K A \ and P±Q = Q, 

det (J - GQK Ai PiG~ l ) = det (J - K Ai QK Ai ) = F(A 2 (x) - x 2 < Vx < t) 

by (2) in [QR12b]. So we are left with estimating 

(/ - QK M Pi)- 1 Qe- aH (K M - I)P 2 1 v -i\ 

(I-P 2 K M P 2 )- 1 P 2 e sH K Ai P 1 J J L2{M 

The last determinant equals det(J — K), with K = Ri y iRi, 2 R 2}2 R 2y i and 

Ri,i = G{I-QK Ai P 1 )- 1 G- x , R 1 , 2 = GQe- sH (K Ai -I)P 2 

R 2 , 2 = (I - P 2 K Ai P 2 )-\ R 2tl = P^KuPtG- 1 . 

Now |det(J -K)- det (7)| < H^Hie 1 ^! 1 < \\K\\ x e 2 for \\it\\i < 1, so the proof will be 
complete once we show that 

(3.11) \\k\\ x < \e- 2 a 2 {t + s )-3 e -f(/3+D 3 / 2 (t+< 

To get this estimate we use (2.1) to write 

II WW <r" II P II II P II II P II II P II 
II-"- 111 S ||-Ki,i||i||-Ki,2||l||-K2,2||l||-K2,l||l 

and use Lemma 3.4, which gives, writing a = 1 + /3, 

\\k\\i < CS~ 3/2 t 1/2 (t 2 + s 2 )3/2 e -|(^ 2 +2< S + S 2 ) 3 / 2 -| f 7 3 / 2 t 3 - S ( ( 7t 2 +2ts+ S 2 )^ 

Now taking s = at we get 

ll^lli < ca^ 2 (l + a 2 f/ 2 t 2 e-^ 1+a ^ 3/2t3 - h ^ t3 , 

where h a (a) = §<r 3//2 + a(a + 2a + a 2 ) - |(ct + 2a + a 2 ) 3 / 2 . Observe that for fixed 
a = (3 + 1 > 4 we have h a (0) = and h' a (0) > 0, which implies that h a (a) > for small 
enough a. Therefore (3.11) holds for small enough a and large enough t, and the result 
follows. □ 

Lemma 3.4. Let R\\, R\ y2 , R 2)2 and i?2,i be defined as in (3.10). Then there is a c > 
such that if t and s are large enough and a = /3 + 1 > 4, 

(3.12a) ||-Ri,i||i < 2, 

(3.12b) ||i? 1)2 ||l < CS-^t-^t 2 + s 2 ) 2 e -2 CTt 3 - S (a < 2 + 2 ts+ s 2 ) ) 

(3.12c) \\R 2 , 2 \h < 2, 

(3.12d) ||i? 2 ,l||l < CS-H 3 / 2 (t 2 + s 2 ) -l/2 e -§( CTi 2 +2is+S 2 )3/ 2 -| CT 3/ 2t 3 +2a4 3^ 

The proof of this result is postponed to Section 4.2. 

4. Estimates of operator norms 

We will use below the following well-known estimates for the Airy function (see (10.4.59- 
60) in [AS64]): 

(4.1) |Ai(x)| < C7e-i x3/2 forx>0, \Ai(x)\ < C for x < 0. 

We start some with some basic integral estimates involving the Airy function: 

Lemma 4.1. There is a c > such that for any m > and a,a',k,t S R such that if 
a > or m > ^a 2 t 2 we have, for large enough t, 

/•rv-i 

-atx A;/„.\2 s- „+2k-l~atm-im 3 / 2 



poo 

/ dx x k e~ atx Ai(x) 2 <ct 2 

J m 
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and 

Ik— 1 „— cdra— |m 3 ' 2 



/*oo /*oo 

/ dx dy x k e~ atx Ai(x + y) 2 e a ' y < ct 2k ~ x e 

Jm JO 



Proof. Using (4.1) the first integral is bounded by 

°° dxx k e- atx -^ /2 = ct 2( - k+ ^ f°° dxx k e-^ + l x3/2 ^. 

The exponent is maximized for x > at x = ^a 2 < m if a < and at otherwise, so 
the first estimate follows from a simple application of Laplace's method, see Lemma 5.1 of 
[CQR12]. The second integral is bounded in the same way after noting that Ai(x + y) < 

_2„3/2_2 3/2 

ce 3 3 s iori,i/>0. □ 

4.1. Estimates used for the upper bound. 

Lemma 4.2. Let Q\ and Q2 be defined as in (2.6). Then there is a c > such that for 
m > —It and t > 1 we have 

||<3i|| 2 <ct 3/4 and ||Q 2 ||2<ct 3/4 . 
Proof. Writing x = 2 1 ' 3 x we have 

/■oo /*oo /'oo /*oo 

\\Qi\\ 2 2 = dx d\e eX Ai(x + m + \) 2 = dxe~ £5: - £rn d\e eX Ai(A) 2 

./0 J ~oa JO J— 00 

/oo 
ciAe A/4 Ai(A) 2 
-00 

by our assumption m > —2t and the facts that e = t _1 and t > 1. Using the estimate 
|Ai(A)| < clA]" 1 / 4 as A -»■ -oo (see (10.4.60) in [AS64]) and (4.1) we deduce that ||Qi||| < 
c£ 3 ' 2 . The bound for Q2 is proved in exactly the same way. □ 

Lemma 4.3. For large enough t > and s > t we have 

dm||Po^, m (2 1/3 -)ll2<e-^ 3+4si 



-it 
and 

-4.3 



dm||P V'-s, m (2 1/3 -)|ll<e-3 J 

2|t| 



Proof. Since |Ai'(x)| satisfies the same bound (4.1) as Ai(x) for x > (see (10.4.61) in 
[AS64]) we have for s > t > 2 and m > —2t that 

||P ^,m(2 1/3 -) Hi = f^ dx 4e 2xs [s Ai(x + m + s 2 ) + Ai'(x + m + s 2 )] 2 

/■oo 

< c(l + s 2 ) / dxe 2:cs -3^+ m + s2 ) 3/2 . 
Jo 

Integrating over m > —2t and scaling m and a? by s 2 we get 

DO /'OO /'OO 

dm ||P ^, m (2 1/3 -)ll2 < c^ 6 / dm / dxe 2- 3 -§(*+ m+ i) 3 / 2 s 3 

2t J-2ts- 2 JO 

= CS 6 / dm dxe 2xs^-l( X+m+ l-2ts^f/^ 

Jo Jo 



OO /"OO 



<cs 6 / dm/ dxe [2,-|( a; + m +i)3/2] s 3 e4tSv ^+^+ T> 
jo Jo 

where in the last line we used the inequality |(x + m + 1 — 2£s~ 2 ) 3 ' 2 > |(x + m + l) 3 ' 2 — 
4is _2 \/x + m + 1 for a? > 0, m > —2t and s > £ > 2. The term in brackets in the first 
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exponential in the last integral is maximized at x = m = 0, and then applying Laplace's 
method as in the proof of Lemma 4.1 gives 

roo 



dm||P VVn(2 1/3 -)ll2<c-'"< ■■' 



making the dependence on m explicit in the notation). Then there is an n > I satisfying 



3„-|s d +4is 
lliii, || j Ufs,m\i ■ 7112 ^ <-• 

/0 

This gives the first bound. The second bound is similar (and slightly simpler). □ 

Lemma 4.4. Let (7™ be the operator defined in (1.6) in [CQR12] for m, L > (here we are 

4 



the following: for fixed L > there is a c > smc/i i/iat /or all m > 



Proof. The proof of this result can be adapted from the proof of Lemmas 5.2 and 5.3 of 
[CQR12]. In that result it is only proved that the above norm is finite, but one can get the 
above estimate by carefully keeping track of the dependence in m. We leave the details to 
the reader. □ 

4.2. Proof of Lemma 3.4. We will use Lemma 4.1 repeatedly without reference. The 
assumption a = j3 + 1 > 4 enters crucially as it ensures in each case that the hypothesis of 
the lemma holds. Throughout the proof we will assume that t > to, where to > 1 should 
be taken as large as needed to make the estimates work. 

Recall the notation introduced in (3.8). In the present case we have a = b = fit 2 , and 
thus recalling that a = I + ft and writing r 2 = at 2 + 2ts + s 2 , M = M ( g t 2 )i and g = Qpt 2 ,t 
to simplify the notation, we have 

p x = p at2 , p 2 = p r2 and Q = P 1 (I + Mg). 

We also define the multiplication operator 

Nf{x) = ifix^fix), 

where, as before, (f(x) = (1 + x 2 ) 1 ' 2 . Finally, we will use repeatedly the decomposition 

K Ai = B P B , 

where, we recall, Bq(x, y) = Ai(x + y). 

Let us start with the first estimate. Since Q = P\ + P\Mg, we have 

(4.2) 
WGQKP^h < WGP^oPohWPoBoPiG^h + \\GP 1 MqB q P N\\ 2 \\N- 1 P B P 1 G-%. 

Now 

/■OO /'OO 

\\GPiBoPoWl = / dx / dy^(xy 2 e- 4tx Ai(x + y) 2 <cr 5 e- 

Jat 2 JO 

while, recalling that Mgf(x) = e 2t ( x ~ at ) f(2at 2 — x), we have 



/'OO /'OO 

\\GPiMqB P N\\1 = dx dy ip(x)- 2 e- 4at3 ki(2at 2 -x + y) 2 ip{y)~ 

Jat 2 JO 

Similarly 



< p-^HAi II 2 \\P,(n- 1 \\ 2 \\c»- 1 \\ 2 < rf- 2 p- A(Tt3 

— e n m iicoH" w ii2ii t / j \\i — cz e 



3„4 C Tt3_4 fr 3/2 t 3 



pco pco 

WPoBoPiG^Wl = dx dyip(y) 2 e Uy Al(x + y) 2 < ct 3 e 

Jo Jat 2 

and one can easily see that the same estimate holds with a possibly larger constant for 
||iV~ 1 i : b-Boi"iG'~ 1 || 2 ,. Putting these estimates together with (4.2) we deduce that 

WGQKPxG^Wx^cte-^' 21 * <\ 
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for large enough t, and then 

Pi.illi < Yl UGQKmPiG- 1 )^! < Y, WGQKa-^G-Yi < 2, 

fc>0 fc>0 

which gives (3.12a). 

We turn now to Ri )2 . Since e~ s (Kai — I) has integral kernel (in x,y) given by 
f_ dX e sX Ai(a; + A) Ai(y + A) (see (2.3) and the paragraph around it), we may use the de- 
composition Qe~ sH (Kai — I)P 2 = QBoPoe s ^BoP2, where e a ^ is the multiplication operator 
defined by (e a ^f)(x) = e ax f(x), so that 

(4.3) Hi^Hi < \\GQB P N\\ 2 \\N- 1 P e^B P 2 \\ 2 . 

Now 

/•oo rO 

\\GPtBoPoN\ll = dx dye- Ux ^x)- 2 Ai(x + y) 2 ip(y)- 2 

J at 2 J — oo 

^ II M l|2 || D ,„-l||2||,„-l||2 „-4crt 3 ^ „ + -2-4at 3 
< ll Al lloolP I** 5 hw h e <cte 

while 

/•oo rO 

\\GP 1 M Q B Q P N\\l= dx dye-^ t3 ^(x)- 2 Ai(2at 2 -x + y) 2 ^(y)- 2 <ct- 2 e- 4at3 

J at 2 J — oo 

in a similar way. On the other hand 

/0 poo 

dx dy<p(x) 2 e 2sx k\{x + y) 2 
-oo J r 2 

ry 
dy e~ 2sy / dx (1 + (x - yf)e 2sx Ai(x) 2 . 



>r* J—oo 

We split the x integral into the regions (— oo, 0] and (0, y\. On the first one we can estimate 
the integral by 

/•oo r0 

||Ai||^/ dye~ 2sy dx(l + (x-y) 2 )e 2sx <cr 4 s- 2 e- 2sr \ 

J r 2 J — oo 

while on the second one we estimate by 

c / dye- 2sy / dx(l + (x-y) 2 )e-l x3/2+2sx < cr A s- l e~ 2sr \ 



Jr 2 JO 

giving 

\\N- 1 P e s tB P 2 f 2 < cr^s^e- 237 - 2 . 
Putting the three bounds in (4.3) gives (3.12b). 

For i?2,2 we observe that ||i-2-^Ai-F2||i < H-P2-B0-P0II2II-P0P0-P2II2, which can easily be seen 
to be bounded by \ f° r large enough t by bounds similar to (and simpler than) those used 
to prove (3.12a), and thus we get (3.12c) in exactly the same way. 

Finally, for R 2t \ we use a similar decomposition as for R\ 2 : using (2.3) we may write 

HPse^AiPiG- 1 !!! < ||P 2 Boe^ /2 Po|| 2 || J Poe^ /2 J Bo J PiG'- 1 || 2 . 
Now 

/"OO /"OO 

\\P 2 B e~ s V 2 Po\\l = / dx dye~ sy Al{x + yf <cs- l r- l e- 

Jr 2 Jo 

and 

/'OO /"OO 

HPoe-^SoPiG- 1 !!^ / dx dy e~ sx Ai(x + y) 2 ip(y) 2 e 4ty <cs-H 3 e 4 

JO Jut 2 

which gives (3.1 2d). 



3' 



-l^^-fa 3 / 2 * 3 
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